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Abstract 

A (d + l)-dimensional dispersionless PDE is said to be integrable if its n- 
component hydrodynamic reductions are locally parametrized by (d — l)n arbi- 
trary functions of one variable. Given a PDE which does not pass the integrability 
test, the method of hydrodynamic reductions allows one to effectively reconstruct 
additional differential constraints which, when added to the equation, make it an 
integrable system in fewer dimensions (if consistent). 
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1 Introduction 



We address the problem of integrability of multi-dimensional dispersionless PDEs of the 
form 

F(u,Ui,Uij) = (1) 

where u is a (vector-) function of d + 1 independent variables. For definiteness, let us 
consider (3 + l)-dimesional PDEs in four independent variables t, x, y, z. Equations of this 
type naturally arise in mechanics, mathematical physics, general relativity and differential 
geometry. Let us look for exact solutions of of the form u = w(i2 x , R n ) where the 
Riemann invariants R 1 , R n solve a triple of commuting diagonal systems 

R\ = \\R)R X) R i y = ^{R)R^ R l z = r ] i (R)R x . (2) 

Notice that the number of Riemann invariants is allowed to be arbitrary! Thus, the 
original mult i- dimensional equation is decoupled into a collection of commuting (1 + 1)- 
dimensional systems in Riemann invariants. In some cases one first needs to rewrite the 
original PDE Q in a quasilinear form to make the method work: see examples in Sect. 2. 
Solutions of this type, known as nonlinear interactions of n planar simple waves, were 
investigated in gas dynamics and magnetohydrodynamics (UEHI- Later, they appeared 
in the context of the dispersionless KP hierarchy [H I fl ^ [TTi l [T^ l 12^ 1 1231 13] . We will call 
a multi-dimensional equation integrable if it possesses 'sufficiently many' n-component 
reductions of the form (j2J) for arbitrary n (the precise definition follows). Exact solutions 
arising within this approach can be viewed as dispersionless analogues of n-gap solutions. 

Partial classification results of (2 + l)-dimensional integrable systems of hydrodynamic 
type (the integrability is understood in the above sense) were obtained in our recent 
publications [121 EH- Integrable second order PDEs of the form u a = f(u xx , u xt , u xy ) were 
classified in [28 . Particularly interesting examples arise in the theory of self-dual and 
Einstein- Weyl spaces |30 | W\ l9| FTT H lllj. in the context of the Dirichlet boundary problem 
in multi-connected domains [22], and the Whitham averaging procedure (in particular, 
the dispersionless limit) applied to (2 + l)-dimensional solitonic PDEs [20J EH EB] • 

We recall, see [32] , that the requirement of commutativity of the flows © is equivalent 
to the following restrictions on their characteristic speeds: 

\i - A* ^ - fi* - rf ' 1 ' 

i 7^ j, dj = d/djtr, no summation! Once these conditions are met, the general solution of 
(j2J) is given by the implicit 'generalized hodograph' formula |32j 

v\R) =x + X i (R) t + /J(R) y + rf(R) z, i = l,...,n, (4) 

where v l (R) are characteristic speeds of the general flow commuting with (j2J), that is, the 
general solution of the linear system 



Substituting u{R}, FT) into ((H) and using (J2J) one arrives at an over-determined system 
for X 1 (R) , fi l (R) , i] 1 (R) and u(R) as functions of the Riemann invariants R l . This system 
implies, in particular, that the characteristic speeds A*,// and t] 1 satisfy an algebraic 
relation which can be interpreted as the dispersion relation for the system (JI}. 

One can show that the requirement of existence of nontrivial 3-component reductions 
is already sufficiently restrictive implying, in particular, the existence of n-component re- 
ductions for arbitrary n ^2]- This phenomenon is similar to the well-known three-soliton 
condition in the Hirota bilinear approach ^HJ (recall that two-soliton solutions exist for 
arbitrary PDEs transformable to Hirota's bilinear form and, therefore, cannot detect the 
integrability) , and the condition of three-dimensional consistency in the classification of 
discrete integrable systems on quad-graphs pQ. One can show that the maximum num- 
ber of n-component reductions a (d+l)-dimensional PDE may possess is parametrized, 
modulo changes of variables R 1 — > f l (R l ), by (d — l)n arbitrary functions of one variable. 
Therefore, we propose the following 

Definition. A (d+1) -dimensional PDE is said to be integrable if its n-component 
reductions are locally parametrized by (d — l)n arbitrary functions of one variable. 

One of the most important examples of PDEs in four dimensions which are integrable 
in this sense is the 'second heavenly' equation, 

Otx + dzy + d xx 6 yy — Q 2 xy = 0, 

descriptive of self-dual Einstein spaces [3U] (see Example 2 of Sect. 2). Its equivalent 'first 
heavenly' form was discussed in the recent publication JT] where it was shown that n- 
component reductions are parametrized by 2n arbitrary functions of a single variable. 
It would be interesting to compare these reductions with the solitonic reductions of [Zj. 
The requirement of existence of n-component reductions parametrized by 2n arbitrary 
functions of one variable appears to be very strong, indeed, the heavenly equation is the 
only nonlinear PDE of the form 

@tx $ zy fi&xxi $xyi $yy)i 

which passes the integrability test (see the Appendix). 

An interesting six- dimensional integrable generalization of the heavenly equation, 

®tl + ®zz + @tx@zy — Qty^zx = 0, 

arises in the context of sdiff(E 2 ) self-dual Yang-Mills equations [SI]. Its n-component 
reductions are parametrized by 4n arbitrary functions of a single variable (Example 5 of 
Sect. 3). We would like to thank M. Dunajski for drawing our attention to this equation. 
Another integrable example is the six-dimensional system 

m t — n x + nm r — mn r , n z = m y + mn s — nm s , 

see Sect. 3. Its n-component reductions depend on 4n arbitrary functions of one variable. 
Under the additional constraints m r = n r = 0, s — x, z = t, this system is descriptive of 
hyperCR Einstein- Weyl structures [Hj; it was investigated in a series of recent publications 

12313 ESI- 
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In the Example 3 of Sect. 2 we apply our method to the four-dimensional PDE 



F = F F + F t 

- 1 xz - 1 - xx- 1 xy ~ - 1 - yt 

which is the second flow of the dispersionless KP (dKP) hierarchy 

F * = -F 2 + F 

F = F F + F + 

1 xz - 1 xx xy ~ - 1 yti 



see, e. g., It is demonstrated that, considered separately, this equation is not integrable 
(as a four- dimensional PDE), the fact which is not at all surprising. What is more 
important, the method of hydrodynamic reductions allows one to effectively reconstruct 
the differential constraint F xt = \F 2 X + F yy which, when added to the equation, generates 
the (2 + l)-dimensional dKP hierarchy. 

Although the method of hydrodynamic reductions provides an infinity of (implicit) so- 
lutions parametrised by arbitrarily many functions of one variable, the question of solving 
the initial value problem for integrable multi-dimensional dispersionless PDEs remains 
open. A detailed investigation of the behavior and singularity structure of solutions de- 
scribing nonlinear interactions of planar simple waves is beyond the scope of this paper. 

2 Examples 

In this section we list some examples of multi-dimensional PDEs which are integrable in 
the sense of hydrodynamic reductions. 

Example 1. Let us consider the first two flows of the dispersionless KP hierarchy, 



F + = -F 2 + F 

1 xt 2 xx * yyi 

F = F F + F + 

^ xz - 1 xx ■ L xy ~ - 1 yti 

which, in the new variables 

F xx = u, F xy = v, F yy = w, F yt = s, Ftf = ^w 2 + w, F xz = uv + s, 

assume the quasilinear form 

U y = V x , Vy = w x , v t = s x , w t = s y , u t = {\u 2 + w) x , u z = (uv + s) x , 

V t = (\u 2 +W) y , V z = (uV + s) y , S x = (\u 2 + W) y , ( \ U 2 + W ) z = ( UV + S ) t ■ 

Notice that u satisfies the dispersionless KP equation, 

(U t - UU X ) X = Uyy. 



(6) 



(7) 
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Looking for reductions in the form u = u(R l , R n ), v = v(R 1 , R n ), w = w^R 1 , R n ), 
s = s(R}, ...,R n ) where the Riemann invariants R 1 satisfy (j2J), and substituting into ((Zj), 
one readily obtains 

diV = [/diii, diS = X l ^diU, diW = (//) 2 9jit, (8) 
along with the dispersion relations 

X = u + (//) 2 , rf = v + 2ufi l + (//) 3 . (9) 
The compatibility condition didjV = djdiV implies 

d^u = -^-^u + -P-1 9 ^ ( 10 ) 
while the commutativity condition Q results in 

dm x = — ^ — :. (11) 
The substitution of (fTTj) into (fTU|) implies the Gibbons-Tsarev system for u{R) and //(-R), 

V"-^. S,9,u = 2^f%, (12) 

i j£ j, which was first derived in [TH| ITH] in the context of hydrodynamic reductions of 
Benney's moment equations. For any solution u of the system (|T2")) one can reconstruct 
A*, rf and v, s, w by virtue of the relations and (fTU|) . which are automatically consistent. 
The general solution of the system (fl"2|) depends, modulo reparametrizations R % — *> f l (R l ), 
on n arbitrary functions of one variable, thus manifesting the fact that PDEs (JBJ) constitute 
a (2+l)-dimensional integrable system. 

Example 2. The so-called second heavenly equation, 

Otx + 8 zy + xx 9 yy — 2 xy = 0, (13) 

is descriptive of self-dual Ricci-flat metrics [3U]. Introducing the variables 9 XX = u, 6 xy = v, 
6yy = w, 9 tx = p, 9 zy = v 2 — uw — p, one can rewrite ()13|) in a quasilinear form, 

u y — v x , u t = p x , v y = w x , v t = p y , 

(14) 

v z = (v 2 -uw - p) x , w z = (v 2 -uw - p) y . 

Hydrodynamic reductions are sought in the form u = u{R}, R n ), v = v (R 1 , R n ), 
w = u^i? 1 , R n ), p = p(R l ,...,R n ) where the Riemann invariants i? 1 ,...,^™ solve a 
triple of commuting hydrodynamic type systems ((21). The substitution into ()14j) implies 

dip = \ l diU, diV = tfdiU, diW = (//) 2 <9iW, (15) 
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along with the dispersion relation 



X 1 = 2vfx l -w- u(ti*y - ^rf. (16) 

Substituting A* into the commutativity conditions ©, and taking into account that the 
compatibility conditions for the relations dip = X l diU, d^v = ^ l diU imply 

didjU = rdiU H — :djU, 

one arrives at the following system: 

r\ j (/Jp—fJ, 1 ) 2 r\ f) i ~ A** )(^7 J —V 1 ) P) 



;i7) 



didjU = 2—, — f : f ■ — rr diudjU. 

1 J rfl — r) l +u(fA3— /j, 1 ) 1 J 



Solving equations (fT7|) for n\rf and u, determining A 8 from (fTfijl and calculating p,v,w 
from equations ()15|) (which are automatically compatible by virtue of ([17))). one obtains 
the general n-component hydrodynamic reduction of the heavenly equation. Moreover, 
the commutativity conditions will also be satisfied identically. 

We emphasize that the system ()17j) is in involution and its general solution depends on 
3n arbitrary functions of one variable. Indeed, one can arbitrarily prescribe the restrictions 
of //* and if to the incoordinate line. This gives 2n arbitrary functions. Moreover, one can 
arbitrarily prescribe the restriction of u to each of the coordinate lines, which provides 
extra n arbitrary functions. However, since reparametrizations K 1 — > f t {R l ) leave the 
system ()17)) invariant, one concludes that general n-component reductions are locally 
parametrized by 2n arbitrary functions of one variable. This supports the evidence that 
the heavenly equation (|13|) is a true four- dimensional integrable PDE. 

Obviously, the same method applies to other equivalent forms of the heavenly equation 
(|13|). For instance, hydrodynamic reductions of the first heavenly equation 

were investigated in detail in [TT]. Another possibility is to work with the evolutionary 
form of the heavenly equation, 

In both cases one can derive analogues of equations ()17j) which, although involutive, look 
somewhat more complicated. 

The requirement of existence of n-component reductions (parametrized by 2n arbitrary 
functions of one variable) is very strong indeed: as demonstrated in the Appendix, the 
heavenly equation (|T3|) is the only nonlinear PDE of the form 



@tx $ zy f(@xxi $xyi @yy) 0? 

which passes the integrability test. 
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Remark. The heavenly equation (and equivalent forms thereof) belong to the class 
of special Monge-Ampere equations which can be defined as follows. Consider a func- 
tion u(x 1 , ...x k ) and introduce a k x k symmetric matrix U = \v,ij\ of its second partial 
derivatives. A special Monge-Ampere equation with constant coefficients is a PDE of the 
form 

M + Mi + ••• + M k = 

where Mi is a constant-coefficient linear combination of all distinct I x I minors of U, 
< I < k. Here, for instance, M is a constant, M k = det U = Hess u, etc. Equiva- 
lently, this PDE can be obtained by equating to zero a constant-coefficient fc-form in 2k 
variables x l ,Ui. It is an interesting problem to classify integrable PDEs within this class, 
in particular, for k — 4. We emphasize that the case k = 3 is understood completely: 
one can show that, for k = 3, any special Monge-Ampere equation is either linearizable 
by a contact transformation (in this case it is automatically integrable by the method 
of hydrodynamic reductions), or contact equivalent to either of the three nondegenerate 
forms 

Hess u = 1, Hess u = un + U22 + U33, Hess u = uu + U22 — «33- 

We have verified directly that these three PDEs are not integrable by the method of 
hydrodynamic reductions. 

Example 3. Let us consider the second flow of the dispersionless KP hierarchy, 

F = F F + F * (18) 

see Example 1. The question is: should it be regarded as a four-dimensional integrable 
PDE? We will see that the answer to this question is negative, moreover, the method 
of hydrodynamic reductions applied to this equations reconstructs additional differential 
constraints which, when added to (fTHj) . generate the (2 + l)-dimensional dKP hierarchy. 
In the new variables 

F%X ^5 XJJ V, Fyy W , Fyl S, -T J), F XZ UV ~\~ S, 

the equation (jl8j) assumes the quasilinear form 

U y = V x , Vy = w x , v t = s x , w t = s y , u t = p x , u z = (uv + s) x , 

(19) 

V t = Py, V z = (UV + S) y , S x = p y , p z = (UV + S) t . 

Looking for reductions in the form u = uffi 1 , R n ), v = v(R 1 , R n ), w = w(R 1 , R n ), 
s = s(R 1 , R n ), p = p(R 1 , R n ), where the Riemann invariants R l satisfy (J2J), and 
substituting into (fl9*]l. one readily obtains 

div = ^diU, diS = XtfdiU, diiv = (fx l ) 2 diU, d^p = X l d i u 1 (20) 

along with the dispersion relation 

rf = v + V + //A\ (21) 
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The compatibility condition didjV = djd{U implies 

djfx 1 n difj, j 



d&u = —^—d iU + ^ : djU, (22) 

3 \1? - n l fj* - fp 3 ' y ' 

while the commutativity condition Q results in 

djlf = d jU , djX = d jU . (23) 

X> - A* 3 3 fjp - 3 

The substitution of (|2*3^1 into (|2*2*j) implies the over-determined system for fi l (R), X l (R) 
and m(-R), 

(24) 

i j, which is analogous to the system (|17j). There is one crucial difference: the system 
({21} is not in involution. Calculating compatibility conditions dk^djii 1 ) — dj{d^ 1 ) = 0, 
dk{dj\ % ) — dj{dk\ l ) = and dk(djdiu) — dj(dtdiu) = 0, one arrives at extra relations 

((//) 2 (A fc - X j ) + (v k ) 2 (\ j - X) + (^') 2 (A l - A fc )) x 

(25) 

(fi k fi j (X k - X j ) + //'//'(,V - AO + ////(A* - \ k )) = 0. 

Thus, there are two cases to consider. 

Case 1. Equating to zero the first bracket in ([25)1 . 

(//) 2 (A fc - A j ) + {fx k ) 2 {\ j - A 1 ) + - A fc ) = 0, 

one obtains A* = a(fi 1 ) 2 + b. The substitution of this ansatz into (|24|) implies (//) 2 <9.,a + 
<9j& = c^-it, so that a = a, b = u + (3 where a,/3 are arbitrary constants. Thus, A 1 = 
ck(^) 2 + m + (3. Substituting X % into the equation dip = X l diU and using (|2T)|) one obtains, 
upon elementary integration, p = ^u 2 + (3u + aw + 7. Expressed in terms of second 
derivatives of F, this constraint reads 

Fxt = -F 2 XX + (3F XX + aF yy + 7. 

One can show that the constants a, (3, 7 are not essentials and can be reduced to a = 
1, /3 = 7 = 0. Thus, the method of hydrodynamic reductions applied to the PDE (jl8j) 
reconstructs the first flow of the dKP hierarchy, F xt = \F 2 X + F yy . 
Case 2. Equating to zero the second bracket in ([25)1. 

/iV (A fc - X 1 ) + /iV'(A J - A*) + /iV fe (^ - A fc ) = 0, 

one obtains A* = + 6. The substitution of this ansatz into (}2*4"j) implies <9.,a + fi J djU + 
fi l (djb + (9ju) = 0, so that a = —v + a, b = —u + /3 where a, (3 are arbitrary constants. 
Thus, A* = + (3 — u. Substituting A* into the equation diS = X l ^ l diU and using (|2*0|) 
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one has, upon elementary integration, s = —uv + au + (3v + 7. Expressed in terms of 
second derivatives of F, this constraint reads 

FxxF X y ~\~ Fyi aF xx -\- (3F X y -\- 7' 

Again, the constants a, j3, 7 are not essentials and can be reduced to zero. The resulting 
constraint F yt + F xx F xy = characterizes stationary points of the flow (|T%|) . We have 
checked that this constraint is integrable in the sense of hydrodynamic reductions (as a 
(2 + l)-dimensional PDE). 

In any case, we conclude that the equation (fTHj) is not integrable as a four-dimensional 
PDE. This is manifested by the fact that the system (|24jl . which governs hydrodynamic 
reductions, is not in involution. 

Example 4. Let us consider the system 

m t — n xi n z — my + mn x — nm x , (26) 

which, in the limit z — t, has been extensively investigated in [27J El Looking 
for reductions in the form m = m(R l , R n ), n = n(R l , R n ) where the Riemann 
invariants R l satisfy (J2J), one obtains 

din = X l dim, // = A*rf — m\ t + n. 

The commutativity conditions (jSJ imply 

djif = djm, didjm = 0, 

hence, up to reparametrizations R 1 —>■ ^(R 1 ), one has 

m = Y,R\ r / i = f(i? l ) + ^^, 

k k 

where f l (R l ) are arbitrary functions of a single variable. The characteristic speeds X 1 
solve the linear system 

A^'-A i P{ Ri )~ fW 

We refer to [2Zj for the general formula for A* and further discussion of this example in 
the (2 + 1) -dimensional limit z = t. Thus, n-component hydrodynamic reductions of the 
system fl2fij1 are parametrized by 2n arbitrary functions of one variable (n functions f l (R i ) 
plus n functions in the general solution of the linear system for A 4 ). Therefore, the four- 
dimensional system (j2T|) is integrable. Notice that it can be obtained as the condition of 
commutativity of two vector fields, 

[d g - md x - Xd x , d y - nd x - Xd t ] = 0, 

A = const, compare with j^j. Some further multi-dimensional generalizations of this 
example are discussed in Sect. 3. 
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Example 5. The six-dimensional generalization of the heavenly equation, 

Qti + ®zz + dtx^zy — dtydzx = 0, (27) 

has been proposed in [3T]- Introducing the variables 9 tx = a, 6 zy = b, 6 ty = p, 6 ZX = q, 
9zz — r ? &tt = pq — ah — r, one can rewrite f)27j) in a quasilinear form, 

a y — Pxi o-z — Qti bt = Pzt b x = q y , b z = r y , q% = r XJ 

(28) 

Pz = {pq -ah- r) y . 

Hydrodynamic reductions are sought in the form a = a(i? 1 , R n ), b = b(R l , R n ), 
p = p(R l , R n ), q = q{R 1 -, R n ), r = r(i? 1 , R n ), where the Riemann invariants 
R 1 , R n solve the commuting equations 

K X = \\R)R\, R i y = ^ i {R)R i z1 Ri = rf{R) R z , R\ = ft(R) i£, R} = Y(R)K- 
The substitution into (J25j) implies 



jd yi ft 

dip = ft dib, dir = —dib, dtq = —dib, d { a = — —8$, (29) 
fi l fi l \i % 

along with the dispersion relation 

rf = fttfq + A> - ftXb - //a - ftrf. (30) 
Substituting \ % into the commutativity conditions 

djX 1 dj[i % djif dj(3 l djY 
AJ - A* ~~ fjj - //* _ rjj -rf~ ft - ft~ 7 J ' - 7* ' 

and taking into account that the compatibility conditions for the relations dtp = ftdib 
imply 



■j 



d&b = -J^dfi + -^—-O.b. 
3 ft -ft ft- ft J 

one arrives at the following system: 

/3J-/3* AJ-A 4 M^-M* 7J-7* ( }(^-^ i )+b(A*-A^)+7 i -7 J ' J ' 
9 (^-/Lt 8 )+6(A i -AJ)+7'-7^ u i uu j u - 



(31) 



Solving equations (|31|) for ft, A 1 , 7* and b, determining if from the dispersion relation 
(|30|) and calculating p, r, g, a from the equations ()29jl (which are automatically compatible 
by virtue of (|31)l). one obtains the general n-component reduction of the equation (|27jl. 
The commutativity conditions will be satisfied identically. We have checked that the 
system (j3Tj) is in involution and its general solution depends, up to reparametrizations 
R % — > (p % (R l ), on 4n arbitrary functions of one variable. 
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3 Mult i- dimensional linearly degenerate systems of 
hydrodynamic type 

In the recent publication [13], we gave a complete characterization of two-component 
(2 + l)-dimensional integrable systems of hydrodynamic type, 



l) t + A { v, W ){ V w )+B M { V w )=* 



which possess infinitely many hydrodynamic reductions. The integrability conditions 
constitute a complicated overdetermined system of second order PDEs for 2x2 matrices 
A and B. In the particular case when the matrix A is assumed to be linearly degenerate, 




A | |. 132) 

these conditions imply 



B 



f(w) 2 fM 

1± — L — aw jLLJ - 

w—v w—v 



f( w ) fM _ av 2 

v—w v—w 



where / is a cubic polynomial, f(v) = av 3 + j3v 2 + jv + 6, and a,(3,j,5 are arbitrary 
constants. A remarkable property of this example is that any matrix in the linear pencil 
B + \iA is also linearly degenerate (that is, reduces to the diagonal form (}3*2*j) after an 
appropriate change of dependent variables). Explicitly, one has 

B = 5B 1 + 1 B 2 + (3B 3 + aB 4 = 

( 1 1 \ / W V \ / W 2 V 2 \ 

W—V W — V \ I ml W — V W — V \ I fi I W—V W — V I I „, 

1 1 ) +7 1 w v I +P I _jP_ J +01 ^ 

V—W V-W / V V-W V-W / \ v - w v - w J \ v - w V - W 

Let us introduce the (5 + l)-dimensional system 

" ) +a( v ) +bJ v ) +bJ v ) +bA v ) +B t ( V ) =0. 

W J, V W J, V W J, V W ), V W J, \ W )r 

Notice that an arbitrary linear combination of matrices A and B%, B 2 , -B3, B4 is linearly 
degenerate! In the new variables m — v + w, n = vw, this systems reduces to 

tth + n x + nm r — mn r = 0, n t + mn x — nm x + m y + n z + mn s — nm s = 0, 

taking a fully symmetric form 

n^i = n x + nrrif — mrif, = rriy + mn~ s — nrris, (33) 
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after the obvious linear change of independent variables. In the limit § — x, rrif = n? = 
it reduces to the system ()2fij) from Example 4. Notice that the system (jHSJ) arises as the 
condition of commutativity of two vector fields, 

[dz — mdg — Xdx + Xmdf, dy — nd s — Xdi + Xndf] = 0. 

Let us demonstrate that the system possesses 'enough' hydrodynamic reductions 
and, therefore, should be regarded as an integrable systems in 5 + 1 dimensions. Looking 
for reductions in the form m = m(R x , R n ), n = n{R}, R n ) where the Riemann 
invariants R 1 solve five commuting systems 

Ri=X i (R)Ri, R\ = ^{R)Ri, Ri = rj i (R)Ri, Rl = ft(R)Rl, Ri = Y(R)Ri, 
and substituting into (JHHj). we obtain 

(A* - nft)d(m = (1 - mft)din, (rf - mj^din = (// - ny l )dim. 
Setting din = ip l dim, one obtains expressions for A* and fi l in the form 
A* = nft + (1 - mft)ip\ fj = nf + (rf - mY)(p\ 
as well as the consistency condition 

OiOjlTL = — — -OiVTi H ; rOjlTL. 

J (p3 - if 1 if 1 - ip'J J 

Inserting the expressions for X 1 and fi l into the commutativity conditions 

djX 1 <9j/i* djrf djft djY 
AJ - A* ~ jJtP - ~ rf -rf ~ ft - ft ~ ^ - rf ' 

one ends up with the following equations for rf, ft, rf, rf 1 and m: 

M. = d iP = = gVzl! 8 m 

rp—rf pi—P* -yJ-Y (l-ml3i)(ri l -rnY)-(l-rnl3 l )(-qi -myi) u j ulj i 

dpV* _ /3'r ? J-7'+m(7'/3 j -7 j /3 t ) Qm (34) 

(pi—ip 1 (1— mfli ){rf — my 1 ) — (1— m/3')(»7-> — m-yi ) 3 ' \ ' / 



OiOj-m — ( 1 _ m/3 j)( r? i_ mr )_( 1 _ m/3 i)(^_ m7 j) OimOjm. 

It has been verified directly that this system is in involution and its general solution de- 
pends, modulo reparametrizations R l — > f l (R l ), on An arbitrary functions of one variable, 
thus manifesting the integrability of the (5 + l)-dimensional system ([33)1 . 
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4 Appendix 

Here we apply the method of hydro dynamic reductions to the classification of integrable 
PDEs of the form 

@tx "I" @zy = fifixxi @xyi @yy)- (35) 

Introducing new variables 9 XX = u, 9 xy = v, 9 yy = w, 6 tx = p, 9 zy = f(u,v,w) — p, one 
rewrites (J35J) in the quasilinear form 

U y = V X , U t = P X , Vy = W X , V t = Py, 

(36) 

v z = (f(u, v, w) - p) x , w z = (/(«, v, w) - p) y . 

Hydrodynamic reductions are sought in the form u = u^R 1 , R n ), v = v^R 1 , R n ), 
w = u^i? 1 , R n ), p = p(R l ,...,R n ) where the Riemann invariants i? 1 ,...,/?™ solve a 
triple of commuting hydrodynamic type systems (J2J). The substitution into (|36|) implies 

dip = \ l diU, diV = fi l diU, diW = (//) 2 <9jW, (37) 

along with the dispersion relation 

X'^fu + fv^ + fM 2 -^- (38) 

Substituting A* into the commutativity conditions (jHJ), and taking into account that the 
compatibility conditions for the relations ()37jl imply 

didjU = rdiU H — :djU, 

one arrives at the following system: 



(39) 



here 

Sij = fun + (j? + l^ j )fuv + ((^) 2 + (l^ j ) 2 )fuw + ^V/w + 

Compatibility conditions for the system (|3*9*|) are of the form 

d k (djiJ) - djtjdkfi*) = (...) djudkU, d^djrf) - d^dkrf) = (...) djud k u } 

dkidjdiu) -djidkdiu) = (...) djud k u, 

where dots (...) denote complicated rational expressions in fi\fi\fi k and rf,^ 3 , rj k whose 
coefficients are functions of the derivatives of / up to third order. Equating these rational 
expressions to zero one arrives at the following system for /: 

fuu fww fuv fwv fvv ^fuw fvvv 0- 

Up to elementary changes of variables, the general nonlinear solution of this system cor- 
responds to the second heavenly equation (fl3~j). 
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5 Concluding remarks 



We have demonstrated that the requirement of existence of 'sufficiently many' n-component 
reductions can be used as the effective criterion providing the test for integrability of 
multi-dimensional dispersionless PDEs. We believe that using the approach outlined in 
this paper along with the available computer algebra packages (we have used Mathematica 
5.0), one can obtain complete lists of multi-dimensional integrable systems within various 
particularly interesting classes, hyperbolic systems of hydrodynamic type being one of 
them. Partial classification results can be found in ^21EHI2E1- The main problems aris- 
ing here are the complexity of integrability conditions (making difficult their geometric 
analysis), and the volume of symbolic calculations required. 

Acknowledgements 

We would like to thank M. Dunajski and M. Pavlov for stimulating discussions and ref- 
erences. 

References 

[f] V.E. Adler, A.I. Bobenko and Yu.B. Suris, Classification of integrable equations on 
quad-graphs. The consistency approach, Comm. Math. Phys. 233, N3 (2003) 513- 
543. 

[2] B. Banos, On symplectic classification of effective 3-forms and Monge- Ampere equa- 
tions, Diff. Geom. Appl. 19 (2003) 147-166. 

[3] A. Boyarsky, A. Marshakov, O. Ruchayskiy, P. Wiegmann and A. Zabrodin, Asso- 
ciativity equations in dispersionless integrable hierarchies. Phys. Lett. B 515, N3-4 
(2001) 483-492. 

[4] M. Burnat, The method of characteristics and Riemann's invariants for multidimen- 
sional hyperbolic systems. (Russian) Sibirsk. Mat. Z. 11 (1970) 279-309. 

[5] R. Carroll and Y. Kodama, Characterization and solution of the dispersionless Hirota 
equations. Nonlinear physics: theory and experiment (Lecce, 1995), 53-59, World Sci. 
Publishing, River Edge, NJ, 1996. 

[6] M. Dunajski, The nonlinear graviton as an integrable system, DPhil thesis, Oxford, 
1998. 

[7] M. Dunajski, L.J. Mason and N.M.J. Woodhouse, From 2D integrable systems to 
self-dual gravity, J. Phys. A: Math. Gen. 31 (1998) 6019-6028. 

[8] M. Dunajski, L.J. Mason and K.P. Tod, Einstein- Weyl geometry, the dKP equation 
and twistor theory, J. Geom. Phys. 37 (2001) 63-92. 



14 



M. Dunajski, A class of Einstein- Weyl spaces associated to an integrable system of 
hydrodynamic type, arXiv:nlin. SI/0311024 

E.V. Ferapontov, D. A. Korotkin and V.A. Shramchenko, Boyer-Finley equation and 
systems of hydrodynamic type, Class. Quantum Grav. 19 (2002) L1-L6. 

E.V. Ferapontov and M.V. Pavlov, Hydrodynamic reductions of the heavenly equa- 
tion, Class. Quantum Grav. 20 (2003) 2429-2441. 

E.V. Ferapontov and K.R. Khusnutdinova, On integrability of (2+l)-dimensional 
quasilinear systems, to appear in Comm. Math. Phys. , jnlin. SI/ 0305044 , 



E. V. Ferapontov and K.R. Khusnutdinova, The characterization of two-component 
(2+l)-dimensional integrable systems of hydrodynamic type, [ami . SI/03 10021 , 

J. Gibbons and Y. Kodama, A method for solving the dispersionless KP hierarchy 
and its exact solutions. II. Phys. Lett. A 135, N3 (1989) 167-170. 

J. Gibbons and S. P. Tsarev, Reductions of the Benney equations, Phys. Lett. A 211 
(1996) 19-24. 

J. Gibbons and S. P. Tsarev, Conformal maps and reductions of the Benney equa- 
tions, Phys. Lett. A 258 (1999) 263-271. 

J.D.E. Grant, On self-dual gravity, Phys. Rev. D 47 (1993) 2606-2612. 

F. Guil, M. Manas and L. Martinez Alonso, On the Whitham Hierarchies: Reductions 
and Hodograph Solutions, rnhn. SI/0209051 

J. Hietarinta, Equations that pass Hirota's three-soliton condition and other tests 
of integrability. Nonlinear evolution equations and dynamical systems (Kolymbari, 
1989), 46-50, Res. Rep. Phys., Springer, Berlin, 1990. 

I.M. Krichever, The averaging method for two-dimensional "integrable" equations, 
Funct. Anal. Appl. 22, no. 3 (1988) 200-213. 

I.M. Krichever, Spectral theory of two-dimensional periodic operators and its appli- 
cations, Russian Math. Surveys 44, no. 2 (1989) 145-225. 

I.M. Krichever, A. Marshakov and A. Zabrodin, Integrable structure of the Dirichlet 
boundary problem in multi-connected domains, arXiv:hep-th/0309010 



V.V. Lychagin, V.N. Rubtsov and I.V. Chekalov, A classification of Monge-Ampre 
equations, Ann. Sci. cole Norm. Sup. (4) 26 (1993) 281-308. 

M. Manas, L. Martinez Alonso and E. Medina, Reductions and hodograph solutions 
of the dispersionless KP hierarchy, J. Phys. A: Math. Gen. 35 (2002) 401-417. 

M. Manas and L. Martinez Alonso, A hodograph transformation which applies to 



the heavenly equation, nlin. SI/0209050 



15 



L. Martnez Alonso, A.B. Shabat, Towards a theory of differential constraints of a 
hydrodynamic hierarchy, J. Nonlinear Math. Phys. 10, N2 (2003) 229-242. 

M.V. Pavlov, Integrable hydrodynamic chains, J. Math. Phys. 44 , N9 (2003) 4134; 
|nlin.SI/03010T0"| 

M.V. Pavlov, Classification of the integrable Egorov hydrodynamic chains, to appear 
in Theor. and Math. Phys. (2003). 

Z. Peradzyhski, Nonlinear plane fc-waves and Riemann invariants. Bull. Acad. Polon. 
Sci. Sr. Sci. Tech. 19 (1971) 625-632. 

J.F. Plebahski, Some solutions of complex Einstein equations, J. Math. Phys. 16 
(1975) 2395-2402. 

J.F. Plebahski and M. Przanowski, The Lagrangian for a self-dual gravitational field 
as a limit of the SDYM Lagrangian, Phys. lett. A 212 (1996) 22-28. 

S.P. Tsarev, Geometry of Hamiltonian systems of hydrodynamic type. Generalized 
hodograph method. Izvestija AN USSR Math. 54, N5 (1990) 1048-1068. 

E.V. Zakharov, Dispersionless limit of integrable systems in 2 + 1 dimensions, in 
Singular Limits of Dispersive Waves, Ed. N.M. Ercolani et al., Plenum Press, NY, 
(1994) 165-174. 



16 



